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Abstract 

We present an exact functional formalism to deal with linear Langevin equa- 
tions with arbitrary memory kernels and driven by any noise structure char- 
acterized through its characteristic functional. No others hypothesis are 
assumed over the noise, neither the fluctuation dissipation theorem. We 
found that the characteristic functional of the linear process can be expressed 
in terms of noise's functional and the Green function of the deterministic 
(memory-like) dissipative dynamics. This object allow us to get a proce- 
dure to calculate all the Kolmogorov hierarchy of the non-Markov process. 
As examples we have characterized through the 1-time probability a noise- 
induced interplay between the dissipative dynamics and the structure of differ- 
ent noises. Conditions that lead to non-Gaussian statistics and distributions 
with long tails are analyzed. The introduction of arbitrary fluctuations in 
fractional Langevin equations have also been pointed out. 



1 



Typeset using REVTpjX 



I. INTRODUCTION 



Noise is a basic ingredient of many types of model in physics, mathematics, economy, 
as well as in engineering. Besides that in each area of research the fluctuations have very 
different origins, in many cases the evolution equation governing the system of interest can 
be approximated by a suitable stochastic differential equation. 

In general the driving forces may be any source of fluctuations and the system can 
be characterized by a given potential. The most popular of those stochastic differential 
equations are the one driven by white Gaussian fluctuations, therefore the problem can 
immediately be reduced to the well known Fokker-Planck dynamics 1-3 . If the fluctuations 
are not Gaussian we are faced with a problem that is hard to solve, but among the different 
types of fluctuations the so call dichotomic noise is a good candidate to study, because in 
general for any potential some conclusions can be drawn 4-8 . If the fluctuations are different 
from the Gaussian ones or the dichotomic noise, in general it is not possible to solve the 
problem for any potential. 

In many situations the technical complications of the model can be diminished by study- 
ing the system in a linear approximation, i.e., around the fixed points of the dissipative 
dynamics. This approximation leads to the study of linear stochastic differential equations 
with arbitrary noises. Besides the simplicity of this kind of equations they were the subject 
of numerous theoretical investigation 9-25 , and also they provide non-trivial models for the 
study of many different mechanism of relaxation in physics, biology and another research 
areas. 

From the previous considerations, it is clear that the usefulness of a linear Langevin equa- 
tion arises from the possibility of working with different kind of noise structures. Therefore, 
one is faced with the characterization, in general, of non-Markov processes. These pro- 
cesses can only be completely characterized through the whole Kolmogorov hierarchy of 
the stochastic process 2 , i.e.: the knowledge of any n— joint probability, or equivalently any 
n— time moment and/or cumulant. 
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In Ref. 26,27 , by using a functional technique, we have been able to characterize arbitrary 
linear Langevin equations with local dissipation, giving therefore a procedure to calculate the 
whole Kolmogorov hierarchy of the process and any n— time moment. Using this previous 
experience, in this paper we are going to generalize our functional technique to tackle the 
more general situation where the dissipative term is non-local in time and the noise is also 
arbitrary; this is what we call a Generalized Langevin Equation, i.e., a linear Langevin 
equation with arbitrary memory and driven by any noise structure 



where £(£), the fluctuation term (i.e., the external noise) is characterized by its associated 
functional, and $ (t) is an arbitrary memory kernel. 

This type of generalized Langevin equation arises quite natural (considering Gaus- 
sian fluctuations) in the context of the Zwanzig-Mori projector operator technique 28,29 ; in 
this case the fluctuation-dissipation theorem is required 30 ' 31 , which imposes that $ (t) = 
(£(t + t)£(t)) jkT . Therefore the dissipative memory must be consistent with the structure 
of the correlation of the Gaussian fluctuations 30 . Nevertheless, if the system is far away from 
equilibrium, the fluctuation-dissipation theorem is not fulfilled and in general the Gaussian 
assumption is not a good candidate to describe the fluctuations of the system. Therefore, 
in a general situation, both the kernel and the noise properties must be considered inde- 
pendent elements whose interplay will determine the full stochastic dynamics of the process 
u(t). Along the paper we will be interested in characterizing this noise-induced interplay 
by assuming different kind of noises and memory kernels, both in the transient as in the 
long-time regime. 

The paper is organized as follows. In section II, after a short review of the functional 
method, we obtain the characteristic functional for a vectorial Langevin equation with mem- 
ory. In section III we apply this result to different situations. First we analyze the case of 
stable noises; after this step, we analyze the stochastic dynamics induced by an exponen- 
tial memory kernel in a process driven by two different fluctuating structures: radioactive 




(1.1) 
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and Poisson noises. Then the interplay between the kernel and the noise properties is em- 
phasized. The necessary noise's properties that guarantee a stationary distribution with a 
long-tail are presented. As an example we introduce the Abel noise that induces this asymp- 
totic behavior. The characterization of a fractional Langevin equation with arbitrary noise 
is also presented. The interplay between the fractional property of the differential operator 
and the Levy noise is analyzed. In section IV we give the conclusions. 

II. FUNCTIONAL CHARACTERIZATION OF ARBITRARY LINEAR 

LANGEVIN EQUATIONS 

A. The characteristic functional method 

In our previous papers, we have presented a complete characterization of Langevin equa- 
tions with local dissipation by means of the characteristic functional of the stochastic process 
u(t) [t G (0, oo)] 

G u ([k(t)}) = (expif~dtk(t)u(t)). (2.1) 

Here k(t) is a test function, and (• • •) means an average over the stochastic realizations of 
u(t) belonging to a given support. The knowledge of the characteristic functional allows a 
full characterization of the process. In fact, from this functional it is possible to calculate 
the whole Kolmogorov hierarchy, and of course any n— time moment and/or cumulant. This 
follows by defining the n— characteristic function G^^kj, tj}™ =1 ) as 

G { :\{k^T 3=1 ) = Gu{[h(t)]), (2-2) 

where the test function kg(t) must be taken as 

k s (t) = h5(t -t 1 ) + ... + kj(t - t n ). (2.3) 

Thus, any n— time joint probability distribution P({uj, t/}" = i) = P n (ui, ti, u 2 , t 2 • • • ; u n , t n ) 
can be obtained by Fourier inversion of G^'\{kj,tj}'j =1 ) 
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PdujMU) = (2^f j dh--- jdk n eM-if^k^G^ak^tj}]^). (2.4) 



3= 

On the other hand, any n— time moment can be calculated as 

d n G^\{u^ =1 ) 



(u(t 1 )u(t 2 )---u(t n )) = {-i) 1 



(2.5) 



dk\dk 2 • • • dk n 

An equivalent formulae holds for the cumulants (or correlation functions), in which case the 
differentiation must be taken from the logarithm of the n— characteristic function 

(2.6) 



{{u(t 1 )u(t 2 )---u(t n ))) = (-i) 



dkidk 2 ■ ■ ■ dk n 

We have showed 26 ' 27 that in order to obtain the characteristic functional of u(t), it is 
only necessary to know the characteristic functional of the noise £ (t) 

G i ([k(t)]) = (e X piJ^dtk(tm). (2.7) 

No any other hypothesis were assumed over the noise. This general formalism allows us to 
deal with arbitrary non-Markovian evolution processes, where there is not a clear underlying 
Fokker-Planck dynamics. 

In what follows we will generalize our functional characterization for the case of a 
Langevin dynamics that includes memory effects, Eq. (1.1). We remark that not any con- 
dition as thermal equilibrium, Gaussian noises, or any other, is imposed on the stochastic 
evolution of the processes u(t). 

B. The multivariable linear case with dissipative memory 

A general case that covers many models of evolution, is the vectorial linear stochastic 
equation 

ju{t) = - £ dt> * (f - f) ■ u(0 + m- (2-8) 

Here, u(t) represent a d— dimensional stochastic process {ui(t)}f =1 , with Ui E (— oo,+oo). 
The memory kernel $ (t) is an d x d matrix with arbitrary kernel functions. The 
d— dimensional vector noise is characterized through its characteristic functional 



G € ([k(0]) = <exp</ oo dtk(0 (2-9) 
where the dot {•} denotes a scalar vectorial product, and k(t) is a vector of test functions 

We want to found an exact expression for the functional of the process u(t), defined as 

G u ([k(t)]) = (expi/ °°dik(i) -u(0>. (2.10) 

The basic idea consists in writing the solution of Eq. (2.8) in terms of the corresponding 
Green function. By denoting the Laplace transform as f(s) = J °° e~ st f(t) dt, from Eq. (2.8) 
we get 

su(s) - u(0) = -<fr (s) • u(s) + £{s). (2.11) 

This allows to express the formal solution for each realization of the noises as 

u(f) = (u(f)) + f dt' A (t - if) ■ m, (2.12) 
Jo 

where we have defined the vector 

(u(*)) o = A(0-u(0). (2.13) 

The d x d matrix A (t) is defined through its Laplace transform 

A(s) = T -L— (2.14) 
si + $ (s) 

where I is a d— dimensional identity matrix and $ (s) is the Laplace transform of the matrix 
kernel $ (t). Eq. (2.14) is equivalent to the evolution 

-^A(t) = - / cft'A(t') •*(*-*'), A(0) = 1, (2.15) 
dt J 

From this equation it is possible to identify the matrix A(t) with the Green function of 
Eq. (2.8). Notice that the vector (u(t)) corresponds to the average value of the process 
u(t) for the case in which the average of the noise is null, i.e., (£(£)) = 0. 
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After introducing the solution Eq. (2.12) into Eq. (2.10), and reordering the time integrals 
Jo ^ Jo dt" — Jq dt" J* dt'Q(t' — t"), where O(t) is the step function, we arrive to a closed 
expression for the characteristic functional of u(t) in terms of the functional of the vector 
noise. Thus, we get 

G u ([k(t)]) = G (u)o ([k(t)])G 5 ([z(t)]), (2.16) 

where we have defined 

G {u)o ([k(t)}) = exp{* f °° dtk(t) ■ (u(t)) }, (2.17) 
and the vector test function z(t) as 

/CO 
dt'k(t')-A(t'-t). (2.18) 

The expression Eq. (2.16) gives us the desired exact functional G u ([k(i)]) as the product 
of two functional. The first one corresponds to the deterministic evolution, or equivalently 
to the averaged process when the noises have null averages. On the other hand, the second 
term comes from the noise characteristic functional, i.e., the stochastic evolution. 

Finally, the n— time characteristic function G^\{kj, tj}™ =1 ) can immediately be evalu- 
ated from the characteristic functional of the noise. The unidimensional case, Eq. (2.2), is 
easily generalized to the vectorial case as 

G^\{k v t 3 }] =1 ) = G u ([k 5 (t)]), (2-19) 
where the vectorial test function kg(t) is 

k 5 (t) = ki<?(f + + k n 5(t - t n ). (2.20) 
Thus, using these last two equations and Eq. (2.16) to (2.18), we get 

G { :\{k v t } n 3=1 ) = exp^EJUk, • (u(f,-)) }G € ([y(f)]), (2.21) 
where the function y (t) reads 
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(2.22) 



Eq. (2.21) gives, in a simple way, the n— characteristic function of the process u(t) by 
evaluating the characteristic functional of the noise with the function y (t). This last function 
is defined in terms of the Green function A(t) of the problem. At this point, it is important 
to remark that our formalism is valid independently of the assumed form for the dissipative 
Green function evolution. Thus, this function may present a monotonous or oscillating 
decay (taking positive and negative values), and in general any one consistent with the 
corresponding kernel. In Appendix A we show a particular 2— dimensional stochastic process 
{ui{t)} 2 i=1 where each component Ui(t) represents the position and the velocity of a memory 
damped harmonic oscillator in presence of an arbitrary noise structure. 



Here we are going to analyze, and characterize with our functional method, many amazing 
situations that arise by choosing different noise structures and a non-local dissipative kernel. 



From now on, we will restrict to a unidimensional case. By assuming a 5— Dirac correlated 
kernel 



III. EXAMPLES 



A. Local Dissipation 



$(0 



lS(t), 



(3.1) 



we arrive to the evolution 



di 



u(t) 



(3.2) 



Therefore, it is possible to obtain A (s) 



l/(s + 7 ). Thus, 



A(t) 



exp [-7 1] . 



(3.3) 



With this solution, it is possible to recapture all the results obtained in Ref. 
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B. Stable noises 



Stable probability distributions play an important role in the theory of sums of random 
variables. This fact follows from the generalized central limit theorem valid for Levy distri- 
butions and the asymptotic power-law decay 32,33 . Here we are going to apply our functional 
formalism with the use of stable noises. 

1. Gaussian noise 

A zero-mean Gaussian noise £(t), with an arbitrary correlation function o"^(r 2 ,ri) = 
(£(t 2 )£(ti)) , is characterized by the functional 2 

(1 poo poc \ 

~2j dT2 J dr iMT2)^(r 2 ,n)A;(ri)J. (3.4) 

Therefore, from Eqs. (2.16) to (2.18) the characteristic functional of the process u(t) results 

(i poo poo \ 

--J dr 2 j^ dnk^k^cr^n)), (3.5) 

where 

<Tu(t2,Ti)=/ dn rfr a A(r 2 - r 5 )^(r 6 ,r a )A(r 1 - r a ), (3.6) 
Jo Jo 

is the correlation function of the Gaussian process u(t). In fact, the n— time characteristic 
function of u(t) is a multivariate Gaussian process. 

2. Levy noise 

Another example that share similar conclusions than using the Gaussian noise, is the 
Levy noise. The characteristic functional of this noise reads 21 

G,([k(t)}) = exp (-y /"^IM^r) , (3.7) 

where < v < 2. Thus, for example, the 1— time characteristic function of the process u(t) 
reads 



G^(k,t) = exp{ik (u(t)} Q }exp 




(3.8) 



where 



E u (t) = a, 



V 



I 



dt A(t) 



V 



(3.9) 



Therefore, memory-like linear Langevin equations driven by stable noises give rise to new 
stable stochastic processes whose correlations properties are defined in terms of the Green 
function A(t) of the dissipative dynamics. Note that this result is valid independently of the 
particular form of the Green function. 



In general, around a fixed point, the linearized dynamics of a dissipative system- 
depending on the parameters of the problem — may presents a relaxation behavior with 
quite different characteristics: let us say, a monotonous or a time-oscillating decay. Here we 
are going to study the relaxation and the fluctuations of the process u(t) by analyzing the 
interplay between these mentioned dissipative behaviors and different noise structures. 

The different behaviors of the dissipative dynamics can be modeled through a non-local 
exponential memory friction 



therefore $ (s) = 5/ {s + A). From Eq. (2.14) we get A(s) = (s + A)/[s(s + A) + 5], which 
after inversion gives 



C. Interplay between non-local friction and noise structure 



$(t) = <Sexp [-At] , 



(3.10) 




(3.11) 



where the frequency wq is 




(3.12) 
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As desired, the dissipative Green function, Eq.(3.11), has a regime where it decays oscillating 
in time. But on the other hand there is a possibility, for 5 < (|) 2 where A (t) decays in 
a monotonous way. Note that the local case (non- memory dissipation), Eq. (3.1), is re- 
obtained considering the limit: A — > oo and 5 — > oo, with 8/X — > 7 (finite). 

Now we are going to present some results using the Green function Eq. (3.11) in con- 
nection with two different structures for the fluctuating term, i.e., a radioactive noise and 
a Poisson white-noise. We will show that the interplay between the deterministic Green 
function and the structure of the noise play a crucial role to determine the relaxation and 
in general any statistical dynamical property of the stochastic process u(t). On the other 
hand, a memory kernel with a power-law decay will be studied in the context of fractional 
derivatives in the sub-section E. 



1. Radioactive decay noise 

The radioactive noise is non- white and Markovian. Its 1— time and conditional proba- 
bility evolve controlled by the master equation 

dP^t) 



dt 



(3 t(e+i)n + iW-enw] 



(3.13) 



Thus, the realizations of this noise start with some initial value £ and at random times 
decrease by finite unit steps until the process reaches the zero value. The constant (3 de- 
fines the probability per unit of time for such discrete steps (f3— decay). The characteristic 
functional reads 2 

& 



(3 J dr exp -p T + i J dr'k(r') 



(3.14) 



This noise clearly does not reach a stationary regime. Using this noise, the 1— time charac- 
teristic function of the process u(t) results 

rt 1 \ & 



GW(M) = e ifc(tt(t)> ° e-^exp 



ik [ rfr'A(r') +0 f dre-^exp ik [ dr'A(r') ) 

Jo J Jo L Jt-T J / 



(3.15) 
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This noise leads to very striking dynamical behaviors for the driven process u(t). In 
the following figures we have plotted, at different times, the corresponding exact probability 
distribution P(u, t). In order to get this object we have anti-transformed Eq. (3.15) by using 
a fast Fourier algorithm. In Fig. 1, the parameters for the memory kernel were chosen as 
5 = 0.1, A = 1. Thus, the Green function decays without oscillations (see inset). The noise 
decay rate is (3 — 1 and its initial value was chosen to be £ = 2. At time t = 1 (full 
line), we see that the probability distribution is highly irregular and also has a 5— Dirac 
contribution (indicated symbolically with the vertical line). This singular behavior arises 
from the particular realization £(t) in which the noise has not yet suffer any decay. 

Due to the non-vanishing value of the noise average, (£(£)) ^ 0, the 5— Dirac contribution 
has different locations at different times. This term can be directly read from the characteris- 
tic functional of the process u(t), Eq. (3.15), which gives: exp (— £o/3t) 5[u — £ Jq dr' A(r')]. 
The discontinuity of the probability distribution P(u,t) is due to the strong change that 
arises in the driving noise after a (3— decay. In fact, in this example, the noise intensity de- 
creases to half of its initial value. At a later time, t — 5 (dotted line), due to the dissipative 
dynamics the probability distribution P(u, t) loss it discontinuous character and seems to be 
accumulated near of the origin. At even later time, t — 15 (dashed line) this accumulation 
seems to increase. As a matter of fact, the stationary distribution is a 5— Dirac centered in 
u — 0. This behavior is a consequence of the coupling between the noise properties and the 
deterministic dynamics of u(t): as the noise intensity dead at long times, all realizations u(t) 
are attracted by the stable point of the dissipative dynamics, which corresponds to u — 0. 

In Fig. 2 we show the case in which the Green function is oscillatory in time (see inset). 
The parameters are 5 = 0.1, A = 0.2, and for the noise we choose (3=1 and £ = 10. Here, 
in contrast to the previous figure, we have increased the initial intensity of the noise, which 
imply that the 5— Dirac contribution decay in a very fast way. On the other hand, this 
higher initial value imply that the "discrete" decay of the corresponding noise realizations 
do not change appreciably the noise intensity. In consequence, at all time the distribution 
P(u,t) results smooth. As in the previous example, all realization of the process u(t) are 
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attracted by the stable point u — 0. Nevertheless, here the transient behavior reflects the 
oscillatory behavior of the Green function. In fact, we note that at successive times, the 
center of mass of the distribution P(u,t) approximately follows the dissipative dynamics of 
the Green function oscillating around the origin u — (see times: t — 1 continuous line, 
t — 5 dotted line, t — 15 dashed line, t = 25 dot-dashed line). In addition, we note that the 
width of the distribution grows approximately up to t = 5. In fact, at this time most of the 
noise realizations have decayed to a null value. For later times, the distribution is mainly 
governed by the dissipative dynamics of u(t). Therefore, its width decreases up to reaching 
a 5— Dirac form in the stationary state. 

We remark that for this noise the stationary distribution P(u) is always determinate by 
the fixed point of the dissipative dynamics. This characteristic is due to the vanishing noise's 
amplitude at the long time regime. In the next example we show a case where the noise 
properties participate, in a crucial way, both in the transient as in the stationary properties 
of the process u(t). 



2. Shot noise 



The stochastic realizations of the shot noise are defined by a sequence of pulses ip(t), each 
one arriving at random independent times U. The characteristic functional of the noise 2 ' 26 ' 27 
£(£), in the interval t G [0, oo], is 



G € ([*(f)])=expQHdT q (r) 



exp [i / k(t)ip(t-T)dt) -1 



(3.16) 

where q(r) is the density of arriving pulses. Therefore, the complete characterization of the 
process u(t) is given by 



exp i 



G u ([k(t)]) = G {u)o ([k(t)])exp^ drq(r) 

where the function ty(t, r) is defined as 

*(i,r)= / dt' A(t - t')^{t' - t) 
Jo 
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roc 

/ k(t)*(t, 
JO 



t) dt 



(3.17) 



(3.18) 



If ip(t) = AS(t) and the density of "dots" q(r) is uniform g(r) = p the noise reduces to 
the more familiar white shot-noise or Poisson noise whose characteristic functional is 

Gt([k(t)]) = exp rft [exp[i^Jfe(t)] - 1]^ • (3.19) 

By adding two statistically independent white shot-noises, with opposite amplitudes, we get 

Gt([k(t)]) = exp ^2p jf dt [cos (Aifc(t)) - 1]^ . (3.20) 

This functional characterizes a noise with a null average, whose realizations consists in the 
random arriving of 8— Dirac pulses with amplitude ±A Note that in the limit A — > 0, 
p — > oo with A 2 p — D/2, this symmetrical white shot-noise converges to a Gaussian white 
noise, i.e., Eq. (3.4) with o"^(r 2 ,r 1 ) = D8{t 2 — Ti). 

Using Eq. (3.20), the 1— time characteristic function of the process u(t) is given by 

G$\k,t) = exp{ik(u{t)) } exp ^2p^ dr [cos (AfcA(r)) - 1]^ . (3.21) 

As we will show below, this expression leads to a rich variety of possible stochastic dynamical 
behavior for the process u(t); both, in the transient as in the stationary state, the different 
behaviors arise from the competence between the different characteristic time scales of the 
Green function [Eq. (3.11)] and those of the noise. In the monotonous regime, 8 < A 2 /4, 
the decay of A(r) can be characterized by the rate 8/X. For 5 > A 2 /4, when w > A/2, the 
oscillatory behavior of A(r) is characterized by the damping rate A/2 and the frequency w . 
On the other hand, the noise is characterized the rate p and the amplitude A. The analysis 
of the behavior of P(u,t) on these parameters can be simplified by defining the rescaled 
process u (t) = u(t)/A, whose evolution can be written as 
du ( 

—±- r L = -8'J^ dr exp[-(r -r )J u {r )+£(r). (3.22) 

Here, we have defined the dimensionless time r = \t and parameter 8' = 5/ A 2 . Furthermore, 
the dimensionless symmetric Poisson noise £ (t') is only characterized by the dimensionless 
rate p = p/\. Thus, the full dynamical properties of the process u(t) are controlled by the 
parameters 8' and p . 
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In the next figures we have obtained numerically P(u, t) from fast Fourier transform of 
Eq. (3.21). In Fig. 3 the parameters of the kernel are 5 = 0.2, A = 1, and for the noise we 
used A = 1, p = 7.5. For these values, the frequency u>o in Eq. (3.12) is complex. Therefore, 
the Green function A(t) decays monotonously in time (see inset). At short times, t = 0.5 
[Fig. 3(a)], the probability distribution presents a series of peaks, each ones separated by a 
distance A = 1. This property is a direct consequence of the nature of the shot-noise, where 
each arriving 5— Dirac pulse produces a shift of magnitude A in the driven process u{t). In 
this short time regime the effect of the dissipative dynamics is negligible. At later times, 
the dissipative contribution becomes appreciable; its effects is to attract all the realization 
of u(t) to the fixed point u = 0. This action, added to the shift effect of the arriving pulses, 
produces an increasing of the width of all the peaks of the distribution; then any value of u 
becomes (approximately) probable. This effect is clearly seen in Fig. 3(b) [t = 1.65]. Note 
that the distribution retained the "sharpy" structure of peaks. This erasing effect is even 
increased with time; see Fig. 3(c) (time t = 2). At later times, all the signature of the peaks 
have disappeared and the distribution goes asymptotically to a Gaussian distribution; see 
Fig. 3(d) [times t = 2.5 (dotted line) and t = 20 (full line)]. 

In general, an asymptotic Gaussian stationary distribution is always expected when the 
rate of the arriving 5— Dirac pulses is greater than the characteristic decay rate of the Green 
function. When the decay of A(t) is monotonous, this condition can be expressed through 
the inequality p > S' , which is clearly satisfied in the previous figure. When this relation is 
not satisfies, deviations from a Gaussian distribution must arise both in the transient as in 
the stationary distribution. We show this effect in the next figures. 

In Fig. 4, the parameters of the Green function (see inset) are the same than in the 
previous figure, [5 = 0.2, A = 1], and for the noise we have used A = 1, p = 0.1. As 
in the previous case, the behavior of P(u, t) arises from the competence between the shift 
effect of the arriving pulses and the dissipative dynamics of u{t). For example, at short 
times [Fig. 4(a), t = 2.5] the probability shows a set of maxima located at ±A with tails 
pointing in the direction to the origin. These tails are mainly originated by the dissipative 
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dynamics that attracts all realizations to the stable point u — 0. On the other hand, here 
the occupation around u — comes from the shift effect produced by two arriving pulses 
with different signs. In contrast to the previous example, at later times [Fig. 4(b), t — 5] 
the evolution of the distribution P(u, t) is mainly governed by the dissipative dynamics and 
the structure of peaks is retained during all the evolution, Fig. 4(c) t = 10, and Fig. 4(d) 
t = 100, where the stationary state is practically reached. Note that in this stationary state, 
the non-Gaussian characteristics are present in u ~ and ±A. 

We note that during all the transient dynamics, there exist a 5-Dirac contribution (in- 
dicated with the vertical lines) centered at the origin, u — 0. This term is originated by 
the noise realization in which not any pulse has arrived up to time t. Therefore, it is ex- 
ponentially damped with a rate 2p, i.e., it can be written as exp (— 2pt) S(u). At each of 
the chosen times, the weight of the 5— Dirac term are respectively: 0.606, 0.368, 0.135 and 
2 x 10~ 9 . Notice that in the stationary state, the 5— Dirac contribution is always completely 
washed out. We remark that, independently of the values of the parameters, this 5— Dirac 
contribution is always present. Nevertheless, note that in the case of Fig. 3, it is rapidly 
attenuated and its contribution is insignificant in the time scale of that plot. 

In the two previous figures we have assumed a monotonous decreasing Green function. 
In general the interplay between an asymmetric Poisson noise, Eq. (3.19), and an oscillatory 
Green function is similar to that founded in the case of the Radioactive noise. In the case of 
a symmetric noise, Eq. (3.21), due to the symmetry of the problem, the transient dynamics 
is approximately similar to that founded for a monotonous decay. Thus, the more relevant 
aspect to analyze is the corresponding stationary distribution. 

In Fig. 5 we study the properties of the stationary distribution P(u) for different values of 
the parameter of the Green function and the noise. In Fig. 5(a) the parameters corresponding 
to the Green function (see inset) are 5 = 2, A = 1, while for the Poisson noise we have taken 
p = 0.3 (dotted line), p = 0.4 (dashed line) and p = 0.8 (continuous line); in all cases taking 
A = y/0A/p. In consequence, the effective amplitude of the noise pA 2 is constant in the 
three graphs. Due to this selection the tails of the distributions coincide, indicating that all 
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the distributions are asymptotically, for large | u \ , Gaussian. This property will be analyzed 
analytically in the next-subsection. 

As in the case of the non-oscillating Green function, Figs. 3 and 4, here we expect a 
Gaussian stationary distribution when the average waiting time between two arriving pulses 
is shorter than the characteristic relaxation time of the Green function. In this case, the 
oscillatory decay of the Green function is characterized by the rate A/2. Thus, we estimate 
that for p > 1/2 a Gaussian statistics arises. This prediction agree very well with the 
plots shown in Fig. 5(a). Nevertheless, note that this inequality does not depend on the 
frequency of the Green function oscillations. In the next figure, we will check the validity of 
this prediction. 

In Fig. 5(b) we show the stationary distribution P(u) by maintaining the noise parame- 
ters and changing the values of the parameters of the Green function. The noise parameters 
were chosen as p = 0.25, A = 1 and for the Green function A = 1, 5 = 0.15 (full line), 
5 = 0.6 (dashed line), and 5 = 50 (dotted line). In the inset we show the corresponding 
Green function for each values. Consistently with our previous analysis, we confirm that 
independently of the value of the characteristic frequency of the Green function, the station- 
ary distribution P(u) is non Gaussian for p < 1/2. Furthermore, we note that by increasing 
the frequency of the Green function, the fast oscillations of A(t) completely drop-out the 
non-Gaussian peaks located at ±A, remaining only one peak centered at u — 0. On the 
other hand, with the chosen values of the noise parameters, a Gaussian distribution only 
arise for a non-oscillating Green function. 

Finally, we want to remark that the criteria for obtaining a stationary Gaussian distribu- 
tion p > 5' , valid when the decay of the Green function is monotonous, and p > 1/2 valid 
for an oscillatory decay, change their validity in a smooth way around the value S' = 1/4, 
which corresponds to the point in which the Green function modify its characteristic decay. 
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D. Long-tail stationary distributions 



Here we are interested in characterizing the stationary distribution of the stochastic 
process u(t); of great importance is to know whether there will be or not a long-tail in the 
distribution. In order to make a general analysis, here we assume that the structure of the 
noise is such that we can write the 1— time characteristic function in the form 

G<£\h, tO = exp dt f (ki A(t))} , (3.23) 

where / (z) is an arbitrary function. Note that this structure is compatible only with a 
white noise. In fact, Eq. (3.23) follows using the test function z(t) = 0(ti — t) A(ii — t) k± 
in the general expression for the functional of u(t), Eq. (2.16), and assuming G^[k(t)] = 
exp J °° dtf(k(t)). Taking the logarithm in Eq. (3.23), and introducing the Laplace transform 
(here denoted as C s ) we arrive to the general expression 

\nG$\k,t = oo) = lims£jlnG'«OM)l = lim£ s [/ (k A(t))] . (3.24) 

If a Taylor expansion of / (z) exists, and we can commute the Taylor expansion and the 
Laplace operator, we arrive to 

oo 

lnG:*(A;) = ^C n A;"£ s [A(0"L =0 . (3.25) 

n=l 

This formulae characterizes the stationary distribution of the stochastic process u(t), if it 
exists, as a series expansion in the Fourier component k. 

From Eq. (3.25) it is simple to see that even when the stochastic transient can be non- 
Gaussian (depending on the driving noise) in the large asymptotic scaling | u | — > oo, the 
Fourier transform of the stationary distribution goes like 

G*(k) ~ exp (i kA - k 2 B + •••), for k ~ (3.26) 

where Ao,B are constants given in terms of the dissipative memory and the structure of 
the noises; then at large scale the behavior is not anomalous. The analysis with a non-white 
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noise is similar and the general conclusion is not changed. So in order to understand the 
occurrence of anomalies or long-tails in the stationary distribution of the process u(t) we 
have to consider driving noises that break the hypothesis Eq. (3.25). This case is achieved 
by structures like the one from the Levy noise, see Eq. (3.7), or its associated one-side 
power-law noise shown in appendix B. 

In the previous sub-section we have shown that the Poisson noise, Eq. (3.19), can give 
rise, during the transient and in the stationary state, to strongly non-Gaussian distributions 
for the process u(t). Nevertheless, this noise can not give rise to long-tail stationary distri- 
butions. As a matter of fact, the characterization of the stationary state induced by this 
shot noise follows straightforward using Eq. (3.24). From Eq. (3.19) its follows 

]nG*(k) = lim p f dr (exp[iAA;A(r)] - 1) , 

t^oo J 

= — {Ei (iAk) — In (A \ k |) - £} , 

7 

where Ei(x) is the exponential integral function 34 , and 8 is the Euler constant. For simplify- 
ing the analysis we have assumed the local Green function Eq. (3.3). Now it is interesting to 
remark about the non-analytic structure of G^(k), perhaps pre-announcing an anomalous 
behavior in the stationary distribution of the process u(t). Nevertheless, it should be noted 
that in the large scale limit k — > this non- analytic structure cancel out, leading therefore 
to an expansion like Eq. (3.26), which imply the absence of any long-tail distributions. This 
is not what happens if the driving noise has a long-tail structure as the one we are going to 
show below. 



Abel noise 

In close connection with a Poisson noise with a random density of arriving pulses, in 
appendix B we have defined the characteristic functional of Abel noise The functional 
of this noise reads [t G (0, oo)] 
2 



G,([k(t)\) 



,jT (l-e ik ®)dtj (^]j a J o (1 -e ik ®)dtj , (3.27) 
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where (z) is the Basset function and r(/i) the Gamma function 34 . Note that in contrast 
with the Levy noise, this functional has well defined integers moments {C(t) q ) if \i > q. Using 
this noise, the 1— time characteristic function of the process u(t) reads 



G$Hk,t) = 2e ^ ( j }) " (J a J\l - e**V) d^j K, ^ a j\l - e^W) dr j . (3.28) 



From this expression (if A(t) > 0,Vt) it is possible to see that the asymptotic behavior of 
the probability distribution is characterized by a one-side power-law distribution of the Abel 
form 

A{trexp{-A(t)/u) 
P(u,t) ~ — i+M — ' l ar g e positive u, (3.29) 



where 

A[t) = a f A(r) dr > 0. 
In the case of local dissipation, using the Green function Eq. (3.3), we get 



(3.30) 



A(t) = - (1 - e'^) . (3.31) 

From this result it is simple to see that the (one-side stable) directed random walk behaves 
like P(u,t) ~ (at)'* / u 1+ ^, as was expected; see appendix B. 

In presence of a non-local dissipation and in the case 5 > (|) , due to the fact that 
the Green function Eq.(3.11) changes its sign oscillating in time, a power-law behavior is 
obtained in a quite unusual way during the transient, i.e., the long-tail changes its support 
from a positive to a negative domain and so on as the time goes on. On the other hand, in 
the stationary regime we get 

A{t = oo) = ^, (3.32) 
o 

then also a power-law is reached at long times, with its long-tail in the positive domain, i.e., 
for u > 0,P(u,*-> oo) ~ (fY u 
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In Fig. 6 we show this amazing behavior for the noise parameters a — 3, fi — ^, 
Ki/2(x) = ^/n/2x exp(— x). In this case we have used for the Green function (see in- 
set) the parameters values 5 = 0.85 and A = 0.2. At time t = 3.5 [Fig. 6(a)] we can see 
that there is a long-tail in the positive domain, but also a sort of exponential decay appears 
in the negative domain as a result of the interplay between the oscillating Green function 
and the Abel noise. At later time t — 5 [Fig. 6(b)] it is possible to see that the long-tail 
has changed its support, showing also an exponential decay but now in the positive domain. 
At time t = 50 [Fig. 6(c)] the stationary distribution is almost reached with the expected 
long-tail in the positive domain. In these figures we have also fitted (with dashed lines) the 
corresponding long-tails with the asymptotic behavior 

P(u,t) ~ 

where sig(u) represents the sign of u. 

To end this sub-section let us call the attention that by subtracting two statistically 
independent Abel noises, it is also possible to obtain the n— time characteristic function of 
a symmetric process u(t), see appendix B. 



r(-/i) 



\A(t)\> 



\siq(u) ■ u 



l+M' 



(3.33) 



E. Fractional derivative evolutions 

Another example that can be characterized with the present functional technique is the 
case of fractional Langevin equations with dissipation. There exist many different way of 
introducing this kind of equations. Here we will consider the evolution 

D? [u (t)} - Mo * - = -yfu {t) + £ (t) , (3.34) 
1 (L — a) 

with 1 > a > 0. This evolution was proposed to simulate the dynamics of financial markets 
where it was found that non-Gaussian driving noises should be the suitable ones 23,24 . The 
proper interpretation of this fractional stochastic equation is actually an integral equation 

u (t) -u = -r? oA">(*)] + oDr a [f (*)] . (3.35) 
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Here, uq is the initial condition, and oD t a is the Riemann-Liouville integral operator 35 

It has been proved 23 that the solution of this equation for each realization of the noises can 
be written as 

u(t) = (u(f)) + f dt'k{t-t')i{t'), (3.37) 
Jo 

where 

(u(t)) = u E aA [-( V t) a ]. (3.38) 

On the other hand, the Green function A (t) is given by 

A(f) = Q(t)t a - 1 E a>a [-{r f t) a ]. (3.39) 

Note that for this model of dissipative fractional equations the average (u(t)) and the noise 
propagate with different Green functions. The generalized Mittag-Lemer function E a> p is 
defined by 

^W = Ef^ a>0,,3>0. (3.40) 



fc=0 



From Eq. (3.37) it is evident that our functional approach allows us to introduce any kind 
of statistics to drive the fractional stochastic evolution. 



Competition between fractional derivative and Levy noise 

Here we are interested in analyzing the competence between the statistic of Levy noise 
and the dynamical effects introduced by a fractional derivative structure. In order to get a 
simpler analysis, here we will assume zero dissipation, i.e., rj — 0. In this case, the previous 
evolution reduce to 

<«(0>o = «o> A W = r|t^- (3 ' 41) 
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Note that from differentiation of Eq. (3.35), these two previous expressions characterizes the 
evolution 

= Dl-°m (3.42) 

We would like to stress that this problem in the frame of our functional approach can be 
done in a simple way. Using Eq. (3.41) and the functional of a Levy white noise, Eq. (3.7), 
the process u(t) is fully characterized by the function [see Eq. (3.9)] 

This expression is only well defined if u(l — a) < 1. Note the restriction on the parame- 
ters {u, a}. This result contrasts with that of a driving Gaussian white noise, where the 
restriction is a E (|, 1]. Thus, in the case of Levy noise, smaller values of the parameter 
a E (0,1] are allowed at the expenses of diminishing the value of v e (0,2] (i.e., large 
Levy-step excursions). 

Finally, we want to remark that Eq. (3.42) driven by Gaussian fluctuations, can be 
mapped with an effective medium approximation in the context of disordered systems 3 ' 36 . 
Note that from our functional approach, it is easy to calculate any correlation functions and 
more complicated objects than the propagator of the system; this is something that is hard 
to obtain in the context of self-consistent approximations 37 . 



IV. SUMMARY AND CONCLUSIONS 

We have completely characterized generalized linear Langevin equations when the usual 
fluctuation-dissipation theorem does not apply. Our central result is an exact expression 
for the characteristic functional of the process u(t), for a general d— dimensional correlated 
process (in appendix A we present a 2— dimensional example) where the memory and the 
noise are arbitrary. In this way we have been able to give a closed expression to get the whole 
Kolmogorov hierarchy, i.e., to calculate any n— joint probability and any n— time correlation 
function or cumulant. 
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We have applied our formalism to many different noise structures, and in this way we 
have shown the noise-induced interplay between the effects of the dissipative memory and the 
structure of the driving noises. Particular emphasis have been put in analyzing the 1— time 
probability distribution P(u,t) and its stationary state; we have shown that the transient 
toward the stationary state strongly depends on the interplay between the dissipation and the 
noise structure. As example we have analyzed the case of an exponential memory function 
in presence of different stochastic driving forces; for example, radioactive noise and Poisson 
noise. Also we have discussed whether or not, at large scale, a non-Gaussian distribution 
appears and when these distributions have long-tails. In order to get this class of stationary 
distributions we have introduced, for the first time, the Abel noise in connection with the 
occurrence of one-side power-law distributions (see appendix B). We remark that this noise 
structure can be useful when studying models where molecular diffusion is an important 
ingredient to be considered. In this context, the relaxation analysis of a generalized Langevin 
particle considering memory dissipation and injection of energy by microscopic random 
contributions characterized by a symmetric power-law, can be also carried out by using the 
present functional formalism. As a matter of fact, asymptotically when the Abel noise is 
symmetric the conclusions, for fi < 2, are in agreement with the structure of the Levy noise. 

Of particular relevance is the use of our formalism to study fractional Langevin equa- 
tions driven by arbitrary noise structures. We have analyzed the competence between the 
noise statistic and the fractional derivative operator. As examples we have compared Levy 
and Gaussian noises and their interplay with the fractional calculus. We remark that the 
possibility of introducing arbitrary noise statistics in fractional Langevin equations is an 
interesting step forward to broad the possible applications of these equations 23 . 

To end this work let us comment that the present functional approach can also be applied 
to the so called delayed Langevin equations. In that case the associated Green function A(t) 
turns to be the crucial ingredient to study different models of delayed Langevin equations 
driven by arbitrary noises. Results along this line will be presented elsewhere. 
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APPENDIX A: A TWO DIMENSIONAL EXAMPLE 



In what follows we will show a multivariable example that can be fully worked out by 
using our formalism. It corresponds to a particle confined in an harmonic potential and 
driven by an arbitrary noise (fluctuating force). The evolution for the position x(t) and 
the velocity v (t) reads 

j t x(t)=v(t), (Al) 

^v(t) = -n 2 x(t) - f dt' $ (t - f) v(t') + £(t), (A2) 
at Jo 

where Q is the characteristic frequency of the harmonic potential. All statistical information 
of the processes [x(t),v (t)] is encoded in the characteristic functional 

G xv ([k x (t), k v (t)}) = (expijdt (k x (t)x(t) + k v (t)v(t))^ . (A3) 

For each realization of the noise, the solution of Eqs. (Al) and (A2) reads 

x(t) = (x(t)) + f dt' X(t - t'W), (A4) 
Jo 

v(t) = (v(t)) + [ dt'd(t-t')^(t'). (A5) 
Jo 

Here, we have defined the averages 

(x(t)) =x(0)+ f dt' (v(t')) (A6) 
Jo 

(v(t)) =v(0)d (t) - n 2 x(0) X(t). (A7) 
The functions $ (t) and X(t) are defined from their Laplace transforms 

^) = ^, g(8)= _f__ , ,A8, 



which is equivalent to 

d 



dt X(t)=m, (A9) 

Igr(t) = -Q 2 X(t) - f dt' $ (t - t') $(t'), (A10) 
dt J 
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with the initial conditions X(0) = and 5(0) = 1. After inserting Eqs. (A4)-(A5) in the 
definition Eq. (A3) we arrive to the exact expression 

G xv ([k x (t),k v (t)\) = G (:c>o ([^(t)])G' < „ )o ([^(t)])^([z(0]), (All) 
where we have defined 

G {x)o (Mt)}) = exp{iJZ°dtk x (t){x(t)) }, (A12) 
G {v)o ([k v (t)}) = exp {ij °°dtk v (t) (v(t)) } . (A13) 
Furthermore, the scalar function z(t) is defined by 



z(t) 



/oo roc 
dt' X(t' - t)k x (t') + J dt'$(t' -t)k v (t'). (A14) 



Thus, knowing the characteristic functional of the noise allows to obtain the full character- 
istic functional of the bidimensional process [x (t) , v (t)]. 

Free Brownian motion: By taking Q = , the previous case reduces in a trivial way to 
the case of a free Brownian particle. In this case the evolution of the particle position is 
given by 

^a=- J f <tf , (t -o^ +{(«). ( Ai5) 

This second order differential equation also arises when modeling a rigid rotator 26 ' 6 . The 
marginal statistic description of this equation can be obtained from Eq. (All) after taking 
k v (t) = 0. Quantities like (cos[x(t)]) or (sin[x(t)]) follow immediately from the real and 
imaginary part of the characteristic functional, and more complex stochastic objects can 
also be calculated analytically. In general, any order differential stochastic equation can 
be analyzed with our functional formalism after introducing additional stochastic process 
corresponding to the different derivatives of the original process. 



APPENDIX B: ABEL NOISE 



Abel was probably the first to give an application of fractional calculus 38 . He used 
derivatives of arbitrary order to solve the isochrone problem in classical mechanics, and the 
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integral equation he worked out was precisely the one Riemann used to define fractional 
derivatives. In a modern context, the Abel- type of integral equation can be written in the 
form 

e 11 P(0 =o"=ir f (J - yY~ l P{y) dy, 11 > 0, (Bl) 



m Jo 

= a"o£>fP(0, 

where a is an arbitrary constant. It is interesting to note that the particular class of nor- 
malized one-side Levy-type of probabilities 

P(0 = ^r^exp (-a/0 , with a > 0, f > 0, (B2) 

are solution of the fractional differential equation Eq. (Bl). The characteristic function of 
the probability density Eq. (B2) can be calculated taking the Laplace transform, then 

G € (*) = <exp(ifcO>=£. [P(0] s =- ik 
2 



(y/-ikaY (2V-ikaj , (B3) 



where K^z) is the Basset function and r(/i) is the Gamma function. Note that if \x G (0, 1) 
the asymptotic behavior of the Basset function shows the expected divergencies for integers 
moments but in general depending on the value of fi we can get finite moments 

<e>=a 9 ^j^, for »>q. (B4) 

We define the characteristic functional of a stochastic process £(£), closely related to the 
power-law distribution Eq. (B2), in the form 



Gt([k(t)]) = ^ \J a £ (1 - e*W) dtj K, jT(l - e *M) dtj . (B5) 

From this expression all the moments of the noise can be calculated. For example the first 
moments read 

m) = a£ r^T' if " > x ' (B6) 

<£(*i)£(* 2 )> = ^^y^ *d - **) + ^n^p if " > 2 - 
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Note that there is a constant term in the correlation function ((£(ti)£(h))} = (£(*i)£(*2)) — 

«(*>)> («**)>, 

<«(«.){(«.)» = j^fc " fa) + ( (^y) 2 ^ if , > 2. (B7) 

As expected this constant term decreases for large values of /i. From now on we will call 
£(t) the Abel noise in honor to that brilliant mathematician. 

Now in order to have a clear meaning of this noise we introduce an alternative interpre- 
tation. Let a stochastic process n(t) be defined by a directed random walk, then for a given 
value of transition rate p we write the master equation 

^^ = p[P(n-l,t)-P(n,t)}, p>0. (B8) 

Using the Darling-Siegert theorem 39 the characteristic functional of the process n(t) can be 
found by solving the masterly equation 

<^hll =ik{t)nU{n,t) + p[U{n-l,t) -U{n,t)], U(n, 0) = 5(n - n ), (B9) 
at 

where k(t) is any test function. The characteristic functional of the process n(t) follows from 
the limit 



G n ([k(t)]) = hm Vn(n,t). (BIO) 

t— >oo — * 

n=0 

Using the generating function U(Z,t) = Y^=o Z n n(n, t), the solution of the generating 
function of Eq. (B9) is, with the abbreviation /* k(t') dt' = K(t), 

U(Z,t) = {Ze iK W) n ° exp QTp ^Ze iK ^~ iK ^ - l] dt^j , (Bll) 

so using Eq. (BIO) we finally arrive to the functional 

G n ([k{t)}) = e iK(co)no exp (pj™ ^ir^dt' _ x j dt ^j _ (B12) 

From this expression it is simple to see that 

G n ([k(t)])=e iK ^ n °Gs([z(t)}), (B13) 
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where G^([z(t)]) is the functional of the Poisson noise (A = 1) evaluated in the test function 
z(t) = J°° fc(t') dt'. This means that the Markov process n(t) controlled by Eq. (B8), with 
a probability rate p, is equivalent to solve a linear Langevin equation driven by a Poisson 
noise £(£) with a uniform density of arriving pulses p. Thus we can write the associated 
stochastic differential equation 

£ = (B14) 

Now, we will assume that the rate p in Eq. (B8) is a random variable with a distribution 
P{p)- This assumption arises naturally in the context of disordered systems. Therefore, the 
final functional can be obtained from Eq. (B12) as 

poo 

(Gn([k(t)]))p {p) =e tK ^ no / G,([z(t)])P(p)dp. (B15) 

J 

This integral can be done for many different distributions P(p), in particular if we use the 
Abel distribution Eq. (B2), we get 

G 6 ([z(t)]) P( P ) d P=^JT)(\J a f {1 - elZit)) dt ) ^( 2 \Z a i°° (1 " e ^ )) ^)' 

(B16) 

where we have used that 1Z e [/ °°(e lz ^ — 1) dt] < 0. This result ends the interpretation 
of the Abel noise as the noisy term appearing in Eq. (B14), in close connection with the 
solution of a directed random walk model with a random (power-law distributed) transition 
rate p. 

Note that the cumulants of the Abel noise are not self-averaging with respect to the 
cumulants of the Poisson noise and the average over P(p). For example using the Poisson 
functional, see Eq. (3.19) with A = 1 [or Eq. (B12)], it is simple to show that its second 
cumulant is ((£(£i)£(£2))) = pS(h — h)- Then the average of the correlation of the Poisson 
noise reads 

f°° T(u - 1) 

I dpP(p)mt 1 Mh)}) = a^^5(t 1 -t 2 ), (B17) 
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result that is different from the calculation of the correlation of the Abel noise, see Eq. (B7). 

Consider now the subtraction of two Abel noises £ c,b (t) statistical independent: £(£) = 
£ c (t) — £ b (t). In this case, we can write the characteristic functional of the process in 
the form 

Gt([k(t)]) = (e W (if~Ht) [?(t) -e(t)} dty 

= Ge([k(t)])Ge,([-k(t)]), (B18) 

where each functional is given by Eq. (B5). From the properties of the Basset function we 
can write this formula in a compact form, using the Kelvin functions, to handle this in a 
more friendly way. Also from this expression, it is simple to see that the first moment of the 
symmetric Abel noise is null, etc. 
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FIGURES 

FIG. 1. Probability distribution, P(u,t), for a memory-like driven process u(t) with an expo- 
nential dissipative kernel and in presence of a Radioactive noise, as a function of u for three 
different times, t = 1 (full line), t = 5 (dotted line), i = 15 (dashed line), taken in arbitrary units. 
The initial condition was chosen as u(t = 0) = 0. The noise parameters are (3 = 1 and £o = 2. The 
parameters of the Green function are 5 = 0.1 and A = 1; the inset shows its monotonous decaying 
behavior as a function of time. The straight line in P(u,t = 1) indicates the 5— Dirac contribution 
(at short times) in the distribution of the process u(t). Asymptotically P(u,t — > oo) goes to a 
5— Dirac located at u = 0. 

FIG. 2. Probability distribution P(u, t) as in Fig. 1 for four different times, t = 1 (continuous 
line), t = 5 (dotted line), i = 15 (dashed line), t = 25 (dot-dashed line), in arbitrary units. The 
noise parameters are (3 = 1 and £o = 10. The parameters of the Green function are S = 0.1 and 
A = 0.2; the inset shows its oscillatory decaying behavior as a function of time. Here, at large 
times the distribution P(u, t — > oo) also goes asymptotically to a (5— Dirac located at u = 0. 

FIG. 3. Probability distribution, P(u,t), for a memory-like driven process u(t) with an expo- 
nential dissipative kernel and in presence of a symmetric white Poisson noise, as a function 
of u for five different times, Fig. 3(a) t = 0.5, Fig. 3(b) t = 1.65, Fig. 3(c) t = 2, Fig. 3(d) t = 2.5 
(dotted line) and t = 20 (full line), taken in arbitrary units. The initial condition was chosen as 
u(t = 0) = 0. The noise parameters are A = 1, p = 7.5. The parameters of the Green function 
are 5 = 0.2, A = 1; the inset in Fig. 3(d) shows its monotonous decaying behavior as a function of 
time. 

FIG. 4. Probability distribution P(u, t) as in Fig. 3, for four different times, Fig. 4(a) t = 2.5, 
Fig. 4(b) t = 5, Fig. 4(c) t = 10, and Fig. 4(d) t = 100 where the stationary state is practically 
reached. The noise parameters are A = 1 and p = 0.1. The parameters of the Green function are 
as in Fig. 3; the inset in Fig. 4(d) shows its monotonous decaying behavior as a function of time. 
The S— Dirac contributions (see text) are indicated with straight lines. 
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FIG. 5. Stationary probability distribution, P(u), for a memory-like driven process u(t) with 
an exponential dissipative kernel and in presence of a symmetric white Poisson noise, as a function 
of u for different values of the parameters for the Green function and the noise. In Fig. 5(a) 5 = 2, 
A = 1, while for the Poisson noise we take p = 0.3 (dotted line), p = 0.4 (dashed line) and p = 0.8 
(continuous line) and A = y/OA/p. In Fig. 5(b) the noise parameters are the same in each plot, 
p = 0.25, A = 1, but we change the parameters for the Green function, A = 1, 5 = 0.15 (full 
line), 5 = 0.6 (dashed line), and 5 = 50 (dotted line). The insets show the Green functions for the 
corresponding different parameters. 

FIG. 6. Probability distribution, P(u,t), for a memory-like driven process u(t) with an ex- 
ponential dissipative kernel and in presence of the Abel noise, £(t), as a function of u for three 
different times in arbitrary units, Fig. 6(a) t = 3.5, Fig. 6(b) t = 5 , Fig. 6(c) t = 50. The initial 
condition was u(t = 0) = 0. The noise parameters are a = 3 and p = \. The parameters of 
the Green function are S = 0.85 and A = 0.2; the inset shows its oscillatory decaying behavior as 
a function of time. The dashed lines correspond to the long-tail fit (see text). In Fig. 6(c) the 
stationary regime of the distribution of the process u(t) has been reached. 
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